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PLASMA INSTABILITIES AND MAGNETIC FIELD GROWTH IN CLUSTERS OF GALAXIES 



in 
o 
o 

(N 

< 

m 

cn 
> 

m 

o 
in 
o 



6 



X 



A. A. SCHEKOCHIHIN,' S. C. COWLEY,^ R. M. KULSRUD,'' G. W. HAMMETT,^ AND P. SHARMA^ 

February 2, 2008 

ABSTRACT 

We show that under very general conditions, cluster plasmas threaded by weak magnetic fields are subject to 
very fast growing plasma instabilities driven by the anisotropy of the plasma pressure (viscous stress) with respect 
to the local direction of the magnetic field. Such an anisotropy will naturally arise in any weakly magnetized 
plasma that has low collisionality and is subject to sturing. The magnetic field must be sufficiently weak for the 
instabilities to occur, viz., /3 > Re'^^. The instabilities are captured by the extended MHD model with Braginskii 
viscosity. However, their growth rates are proportional to the wavenumber down to the ion gyroscale, so MHD 
equations with Braginskii viscosity are not well posed and a fully kinetic treatment is necessary. The instabilities 
can lead to magnetic fields in clusters being amplified from seed strength of ~ 10"'* G to dynamically important 
strengths of ~ 10 //G on cosmologically trivial time scales 10* yr). The fields produced during the amplification 
stage are at scales much smaller than observed. Predicting the saturated field scale and structure will require a 
kinetic theory of magnetized cluster turbulence. 

Subject headings: galaxies: clusters: general — instabilities — magnetic fields — MHD — plasmas — 
turbulence 



Plasma in clusters of galaxies is likely to be in a turbulent 
state. Cluster mergers produce bulk flows at scales L ^ I Mpc 
with velocities U ~ 10^. ..10^* km s"'. Assuming an intracluster 
medium (ICM) with temperatures T 10^...10* K, these bulk 
velocities are comparable to the thermal (or sound) speed Vth = 
(r/m,)'/^ (here m,- is the ion mass and, for simplicity, we take 
the ICM to be made of hydrogen plasma with equal electron 
and ion temperatures). Flows of this magnitude can drive tur- 
bulence in the ICM. The characteristic time for the turbulence to 
be established is L/U ^ 10^ yr. Numerical simulations of clus- 
ter evolution (Norman & Bryan 1999; Roettinger et al. 1999; 
Ricker & Sarazin 2001) support this qualitative picture. In re- 
cent years there has been a rising interest in measuring the clus- 
ter turbulence both with future instruments {Astro-E2; see Inog- 
amov & Sunyaev 2003) and with present ones {XMM-Newton; 
see Schuecker et al. 2004). 

The ICM is fully ionized. The ion kinematic viscosity is 
V ^ v^/vii, where v-,; = 47rne'*lnAOT,. ^/"^j-^l^ is the ion-ion col- 
lision frequency, n ^ 10"^... 10"^ cm"^ is the ion number den- 
sity (for cluster cores), e is the electron charge, and In A ~ 20 is 
the Coulomb logarithm (e.g., Helander & Sigmar 2002). This 
gives Reynolds numbers in the range Re ^ 10^... 10^. As the 
turbulent velocities at the outer scales are < v'th, the turbu- 
lence in the inertial range will be subsonic. It is natural to 
assume that Kolmogorov's dimensional theory should apply at 
least approximately. Then the viscous scale of the turbulence 
is li, ^ LRe""'''"' ^ 10. ..30 kpc. These numbers appear to agree 
quite well with observations of turbulence in the Coma cluster 
via pressure maps (Schuecker et al. 2004). 

If we estimate the magnetic diffusivity of the ICM using the 

standard Spitzer (1962) formula ?; ^ T~^l^mJ'^ e^c?' In A/47r, we 



MHD turbulence with large Prandtl numbers, we know that un- 
der these conditions a small-scale dynamo operates: magnetic 
fluctuations are amplified by random stretching of the field lines 
(Jaffe 1980; Roland 1981; Ruzmaikin et al. 1989, see a recent 
account of the theory and simulations of small-scale dynamo by 
Schekochihin et al. 2004 and references therein). The magnetic 
energy grows exponentially at the turnover rate of the viscous- 
scale eddies: the exponentiation time is (L/f/)Re~'''^ ^ 10* yr 
The resulting fields have a folded structure: they are organized 
in long thin flux sheets (or ribbons) with direction reversals at 
the resistive scale and field lines remaining relatively straight 
up to the viscous scale or, in the nonlinear regime, even to the 
outer scale of the turbulence. 

This picture is hard to reconcile with the observed magnetic 
fields in clusters. The dynamo-generated fields are expected to 
have reversals on subviscous scales down to the resistive scale, 
which is /,j ~ Rm~'''^L ~ 10"*... 10^ km in clusters, a tiny dis- 
tance. Published rotation measure (RM) data for clusters sug- 
gest tangled fields with B ~ 1...10 /iG and a typical reversal 
scale on the order of 1 kpc (Feretti et al. 1995; Clarke et al. 
2001; Taylor et al. 2002). Recently, Vogt & EnBlin (2003) used 
RM measurements from extended radio sources in several clus- 
ters to extract magnetic-energy spectra. The spectra peak at 
^ 1 kpc and decrease at smaller scales, so one can reasonably 
take 1 kpc to be the resolved characteristic scale of field varia- 
tion. While this is about an order of magnitude below the vis- 
cous scale, it is certainly much larger than the resistive scale. 

These apparent inconsistencies between theory and obser- 
vations are a serious problem. The generation of direction- 
reversing folded fields is a fundamental property of random 
shear. It does not depend on the particular character of the tur- 
bulent flow and does not require very large Re. Even if the clus- 
ter fields originally owe their existence to some external mech- 
anism (e.g., Kronberg et al. 2001) rather than to in situ gen- 
eration by turbulence, any magnetic field introduced into the 



find an extremely small value leading to enormous magnetic 
Reynolds numbers Rm ^ 10^"^... 10^'. The magnetic Prandtl 
number Pm = Rm/Re for the ICM can, therefore, be as large 
as 10^'^. From theory and numerical simulations of isotropic 
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ICM is tangled by turbulence and rendered indistinguishable 
from a small-scale-dynamo-generated field on the timescale 
of ~ 10^ yr. 

The discrepancy between the predictions of MHD models 

and the observed cluster turbulence has prompted us to reex- 
amine the MHD approximation and to recognize that it is not, 
in fact, appropriate for the ICM. For the cluster plasma, the 
ion mean free path greatly exceeds the resistive scale: A^fp ~ 

Vtk/vii ^ Re~'/"*/^ ^ 1..10 kpc. Collisionless effects are impor- 
tant at scales below Amfp, and the fluid MHD description is not 
vahd. This fundamentally alters the dynamics of both velocity 
and magnetic fields at these scales. 

It can be shown that at frequencies below the ion cyclotron 
frequency fi, = eB/cm, and scales above the ion gyroradius p, = 
Vth/fi(, the equations for the fluid velocity u and magnetic fleld 
B have the following general form 

du ( r" T B-VB 



dt 



B Vu-BV u, 



(2) 



where d/dt = d/dt + u-V,p± and p\\ are the perpendicular and 

parallel plasma pressures, respectively, and b = B/B. We have 
dropped the diffusion term in Eq. (2). Equation (1) is valid pro- 
vided the ions are magnetized, i.e., p, <C Amip. For the ICM, 
this requirement is satisfled if B » 10"^^ G, which is far below 
the observed field strengths of 1...10 /iG. The lower limit for 
dynamically important fields is the field strength corresponding 
to the energy of the viscous-scale eddies (~ pv^Re"'''^), which 
gives B 10 fxG. Thus, the ICM would be very weU magne- 
tized already for dynamically weak fields. 

In such a plasma, the fundamental property of charged parti- 
cles moving in a magnetic field is the conservation of the first 
adiabatic invariant iJb=mv\/2B. When X^f^ p„ this conser- 
vation is only weakly broken by collisions. As long as /i is con- 
served, any change in B must be accompanied by a proportional 
change in p±. Thus, the emergence of the pressure anisotropy 
is a natural consequence of the changes in the magnetic field 
strength and vice versa. The following heuristic argument re- 
veals the connection. Summing up the first adiabatic invariants 
of all particles, we get p_l/B = const. Then 



1 dpA 



1 dB 
B~dt 



P±-P\\ 



(3) 



p± dt 

where we have assumed that p = const and the pressure 
anisotropy p±-p\\ is relaxed by collisions at the rate i^u and 
remains small compared with the total pressure p = pv^. From 
Eq. (2), the field strength evolves according to 



1 dB 
B dt 



Vu 



(4) 



(assuming incompressibility V • u = in accordance with the 
earlier observation that the inertial-range motions are subsonic). 
Using this and Eq. (3), we get p±-p\\ ~ (pv^i,/i/,i)bb : Vu. 

A more formal kinetic calculation of the pressure anisotropy 
can be done starting from the kinetic equation for the magne- 
tized ions (to the lowest order in kpt; see Kulsrud 1983) 

V||B\ df_ 
dv± 



dt " " 2 



bb : Vu- V-u + vii 



B 



b 



du „ 

--fV||V|,u 



VhB' 



2 B 



^ = C[f,f], 

dvii 



(5) 



where V|| = b • V, vj^ and V|| are velocities perpendicular and 
parallel to B and the right-hand side is the collision term. A 
perturbation theory is constructed under the assumption that the 
fluid velocity is much smaller than the thermal speed and varies 
at scales much longer than the mean free path. This is appropri- 
ate for the viscous range of Kolmogorov turbulence in cluster- 
type plasmas. Assuming Kolmogorov scaling, the velocity at 
the viscous scale is m ~ Re"'/Sth. The wave number associated 
with the viscous scale is Re^^^L"^ ReT^^^X^. We in- 
troduce the small parameter e ~ A:,^Amfp ~ u/v^ ~ ReT^^'^. In 
cluster plasmas, e ^ 0.1. ..0.3 — while this is not really very 
small, it is convenient to order all terms in Eq. (5) with respect 
to e, which we believe to be the essential physical small pa- 
rameter in the problem. Namely, we assume Amfp ~ e and 
d/dt'^ Vm ~ e^vth/Amfp (this orders out compressible motions, 
so V • u = 0) and take the zeroth-order distribution function 
/o to be a Maxwellian. For simplicity, we use the Lorentz 
pitch-angle-scattering form of the linearized collision operator 
C[/o,/] = (l/2)z/,-,(v,h/v)3(a/90(l -e)df/dt where ^ = v,|/v 
and the derivatives are at constant v = (v^+vp^/^ (e.g.. He- 
lander & Sigmar 2002). Then the anisotropic ion distribution 
function up to order is 



/(V±,V||): 



n exp(-v^/2v^t,) 
(27rv2)3/2 



1 + 



bb: Vu v^(v2-3v 



.(6) 



Computing p±^ and p\\, we obtain the result that is the lowest 
order (in kpi) term in the plasma pressure tensor first derived by 
Braginskii (1965) (who used a somewhat different ordering in 
his perturbation scheme) 



p^-p\\ = 3pvB^h 



Vu = pxiX 



(7) 



where A ~ is the dimensionless measure of the pressure 
anisotropy and is the Braginskii ion viscosity, which is anal- 
ogous to the ordinary viscosity for the unmagnetized case: ^ 
vthAmfp ^ ^th/ '^ii- It is this viscosity that we used in our defini- 
tion of Re and the viscous scale. For simplicity, = 7] = T, so 
the electron contribution to the pressure anisotropy is subdom- 
inant by a factor of (nielm^)^!'^ ~ 0.02. 

Equations (1), (2), and (7), together with V • u = 0, are a 
closed set. The Braginskii viscosity has two key properties. 
First, the velocities that it dissipates are those that change the 
strength of the magnetic field [Eq. (4)], so the small-scale dy- 
namo action is associated only with velocities above the viscous 
scale. Second, velocity gradients transverse to the magnetic 
field (e.g., shear-Alfven-polarized fluctuations) are undamped. 
Velocity fluctuations can now penetrate below the viscous cut- 
off. In what foUows, we will show that the Braginskii viscosity 
not only fails to damp all kinetic energy at the viscous scale but 
triggers fast-growing instabilities at subviscous scales. 

It has, in fact, been known for a long time that anisotropic 
distributions lead to instabilities (Rosenbluth 1956; Chan- 
drasekhar et al. 1958; Parker 1958; Vedenov & Sagdeev 1958). 
A vast literature exists on these instabilities and their (mostly 
space physics) applications, which we do not attempt to review. 
The standard approach is to postulate a bi-Maxwelhan equi- 
Ubrium distribution with 7j_ (e.g., Gary 1993; Ferriere & 
Andre 2002, and references therein). We do not need to adopt 
such a description because Eqs. (1), (2), and (7) incorporate the 
pressure anisotropy in a self-consistent way. 

Let us imagine a turbulent cascade that originates from the 
large-scale driving, extends down to the viscous scale, and gives 
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rise to velocity and magnetic fields u and B that vary on time 
scales > (Vm)"' and on spatial scales > krj . We study the sta- 
bility of such fields. The presence of turbulent shear (velocity 
gradients) gives rise to the pressure anisotropy given by Eq. (7). 
We now look for linear perturbations 5vl, 5B, 5p±, and 5p\\ that 
have frequencies lo ^ Vm and wavenumbers ki,. With re- 
spect to these perturbations, the unperturbed rate-of-strain ten- 
sor Vu can be viewed as constant in space and time."* Lineariz- 
ing Eq. (1) and neglecting temporal and spatial derivatives of 
the unperturbed quantities, we get 



du± oc k_L, and the dispersion relation is 



2kl 
' kl 



LO +iujv2'K 



k\v(h 



-iupSu = -jk5 ( p± + 



+ h5 



SttJ 
V|| {p±-p\ 



+ [P±-P\ 



'-( 



p±-pr 



(5(b- Vb) 



ViiB' 



B 



(8) 



[k\-2k\^i^ + 2k^Pr^ v^ = 0. (11) 

In the weak-field limit (/? | A |"' ), the third term leads to insta- 
bility. likyk\ < |A|i/2 (the case of parallel or nearly parallel 
propagation), the second term is negligible and, for A < 0, we 
get 7 ~ ^||Vth|A|'/^, as in the case of the firehose instability. 

When ki_ is not small (oblique propagation), the second term 
in Eq. (11) dominates the first. It contains the effect of the col- 
lisionless kinetic damping due to resonant wave-particle inter- 
actions (Barnes 1966) (a physical discussion of the role of res- 
onant and nonresonant particles in the mirror instability can be 
found in Southwood & Kivelson 1993). The growth rate is 



where 5 (b • Vb) = ikndh and from the linearized Eq. (2), 5h = 
— (fc||/a')(5u^ and 6b/b = k_L • 6u±/u;. In the resulting disper- 
sion relation, it is always possible to spht off the part that corre- 
sponds to the modes that have shear- Alfven-wave polarization, 
5u oc k_L X b. 



7 = 



1/2 



l^iikth 



A 1- 



'^11 
2kl 



1 + - 



(12) 



w2 = 4v2 (A + 2/ri), 



(9) 



Modes with > \/2fc± are unstable if A < (slow-wave- 
polarized firehose instability); modes with k± > /cy / V2 are un- 
stable if A > (mirror instability). The growth rate in either 
case is 7 ^ fc||Vth|A| ~ e^fc|| Vth. 

The instabilities described above have growth rates propor- 
tional to . Thus, MHD equations with Braginskii viscosity 
do not constitute a well-posed problem. This means, for exam- 
ple, that any numerical simulation of these equations will blow 
up at the grid scale unless a small additional isotropic viscos- 
ity is introduced. The linear-in-A:|| behavior of the growth rate 
is only modified if the finiteness of the gyroradius pi is taken 
into account. This is done using the general plasma disper- 
sion relation (e.g., Davidson 1983) with the equilibrium distri- 
bution function (6). The growth rate of slow- wave-polarized 
instabihties [Eq. (12)] peaks at kpi ~ 1 with 7max ~ Afi,-. For 
the shear-Alfven-polarized firehose, the fastest growth is for 
k±/k\\ ~ y/2/3, and the peak growth rate 7n,ax ^ lAp/^fi,- is 
achieved at kpi ~ 1. We omit the derivation of these results. 

We have shown that, given sufficiently high /3, firehose and 
mirror instabilities occur in the regions of decreasing (A < 0) 
and increasing (A > 0) magnetic field strength, respectively. 
The viscous-scale motions associated with the turbulent cas- 
cade from large scales will stretch the magnetic field and 
thereby produce regions of both types: the typical field struc- 
ture resulting from random stretching is folded flux sheets (or 
ribbons) with field amplification regions containing relatively 
straight direction-alternating fields and curved corners where 
the field is weaker (the field curvature and strength are anti- 
correlated; see Schekochihin et al. 2004). If B > lO"'^ G, the 
plasma is magnetized (p, < A^fp) and this structure is intrin- 
sically unstable: straight growing fields to the mirror, curved 
weakening fields to the firehose instability. Since the instabili- 
ties are much faster than the turbulent stretching, their growth 
and saturation will have a profound effect on the structure of the 
field. We do not yet have a detailed theory of this process and 
only give a qualitative discussion of what can plausibly happen. 

When turbulent stretching acts on some initial weak field Bq, 
pressure anisotropics arise, Aq ~ {l/i'ii)dlnBo/dt [seeEq. (3)], 
and drive mirror and firehose instabihties. Both produce fluc- 
tuating fields SB curved at scale pifi ^3x 10^(10"'** G/Bq) pc. 
Let us concentrate on the mirror instability as it amplifies 

Balbus (2004) recently considered an instability triggered by Braginskii viscosity when the unperturbed Vu consists of rotation and radial shear in a Keplerian disk 
(see also Sharma et al. 2003). His instability had a growth rate 7 ~ | Vu|. In contrast, we are interested in instabilities that are much faster than the rate of strain and 
arise when turbulent stretcliing is present. 



where /3 = 2v\lv\ = Zirpv^jB'^. When the magnetic energy is 
larger than the energy of the viscous-scale eddies, /? <C | A|"' ~ 
Re'^^, Eq. (9) describes shear Alfven waves. In the opposite 
limit, /3 » |A|"\ and if A < 0, an instabiUty appears with the 
growth rate 7 = fcuVthlAI'/^ ^ e^H Vth. Since k\\ ^ k^, the insta- 
bihty is faster than the rate of strain, 7 ^ Vu ~ ek^vib, in ac- 
cordance with the assumption made in the derivation of Eq. (9). 
This instabihty, triggered by p\\ > p±, is called the firehose 
instability. It does not depend on the way the pressure pertur- 
bations are calculated because it arises from the perturbation of 
the field- line curvature b • Vb in Eq. (8): to hnear order, it en- 
tails no perturbation of the field strength and therefore does not 
alter the pressure. 

In order to determine the stability of perturbations with other 
polarizations, Sp± and (5/?|j have to be computed. We retain the 
assumption that lo ;» Vm and consider scales smaller than the 
mean free path, ^A,nfp » 1. At these scales a fully kinetic de- 
scription must be used. It is sufficient to look for "subsonic" 
perturbations such that iv <C A:vth because the high-frequency 
perturbations are subject to strong coUisionless damping and 
cannot be rendered unstable by a small anisotropy. We linearize 
Eq. (5) around the distribution function (6) and calculate dp± 
and from the perturbed distribution function while neglect- 
ing terms of order (w/A:|| Vth)^ and higher. This leads to 



: klvl I 



+ ki 



A+2p-^)6u±- 



b ^ k • J 

fell 



2{-A + (r')vl-iw 



2.^ + 2" 

^,,1 kj 



(10) 



k_L • i5u_L. 



For the modes with shear-Alfven-wave polarization, we re- 
cover Eq. (9). If we now dot Eq. (10) with b, we get (5m|| = 
-k^ • Su±/k\^^. Therefore, the perturbations with (5m|| y^O are in- 
compressible and slow-wave polarized. For these perturbations. 
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field strength in the linear order The amplification rate is 
70 ^ Aofi/.o ^ 10"** (Bo/ 10"'** G) yr"'. Suppose for a moment 
that the instabiUty does not saturate until the field strength in- 
creases to some value Bi ~ Bq + 5B, where 8B/Bq ~ 1. Then 
the field at scale /9, o is unstable to perturbations at scale p,- 1 ^ 
(Bo/Bi)piX) < pix). This is a secondary mirror instability that 
feeds on the pressure anisotropy due to the increasing field Z?i 
excited by the primary. Braginskii theory cannot be used to 
calculate this anisotropy because the field Bi is curved on a col- 
lisionless scale pn). However, p is still conserved at this scale 
(conservation is broken at scale and we can adapt the esti- 
mate (3) to the coUisionless regime by assuming that pressure 
anisotropy is relaxed in the time particles streaming along the 
field line need to cover the distance p, o : we replace f ,, in Eq. (3) 
by Vth/p;.o = and write Ai (l/f^aO-^lnBi/fif ^ lo/^u) ^ 
Aq. Thus, the pressure anisotropics driving the secondary and 
the primary instabilities are the same. The growth rate of the 
perturbations at scale pi^i is then 71 ~ Aof2,-,i ~ {Bi/Bq)^q. 
This argument can be iterated. As the growth rate 7max oc B, 
magnetic field grows explosively, dB/dt ocB^, until it is strong 
enough to cancel the pressure anisotropy and shut down the in- 
stabilities: /3 ~ A-i ~ Re^/^. This gives B ~ 10 ^G, which is 
the lower bound for dynamically important fields and is com- 
parable to the observed field strength in clusters. Thus, a seed 
field of order 10"^^ G (attainable by primordial mechanisms; 
see, e.g., Gnedin et al. 2000) can be amplified to dynamical 
strength in about 10^ yr. 

The above argument depends on the assumption that fluctu- 
ating fields would grow IoSB/Bq^ 1. What if the saturation of 
the instabilities is quasiUnear, with |(5BP/Bq <C 1? In the quasi- 
linear theory, small fluctuations scatter particles and effectively 
increase their collision frequency to i/ql ^ (|(5BP/Bq) 7j„ax- 
The pressure is isotropized and the instabilities are quenched. 
Since the effective viscosity (~ v^/pql) of the plasma is re- 
duced, the Reynolds number increases Rcql ~ (^/vth) i^ql oc B 
and the turbulent cascade extends to smaller scales, giving rise 
to faster field stretching (by the viscous-scale eddies at the rate 
oc Rcq^^) and to pressure anisotropy A Rcql^. This argu- 
ment can be formalized somewhat and gives explosive growth 



1 /2 

of both Rcql and B until f3 ^ Rcql • At this point the insta- 
bilities should start shutting down with Re dropping back to its 
original value based on particle colUsions and magnetic field 
following up so that f3 ~ Re'^^. The characteristic time for this 
process to complete itself is ~ 10^ yr Thus, both the quasilin- 
ear scenario (SB/Bq 1) and the case of JB/Bq ~ 1 appear to 
produce a fast amplification stage with similar end results. 

After this stage, f3 < Re'^^, the field can be stretched with- 
out going unstable and the folded structure could finally be set 
up. However, differences remain between the MHD model with 
isotropic viscosity (Schekochihin et al. 2004) and the ICM. 
First, the Braginskii viscosity implies that a cascade of shear 
Alfven waves may exist between the viscous scale and the ion 
gyroscale. The properties of this cascade remain to be under- 
stood. Second, even when the rms field strength Bmis is suf- 
ficient to quench the instabilities, there will always be regions 
where the field is weak. For the folded fields, these are the 
bending regions (curved field). They can become unstable un- 
less the field-line curvature there is larger than 1 / p,,bend> where 
A.bend is the gyroradius computed with the local field strength 
Bbend ^ Brms and, therefore, is larger than p, based on Br^s by 
a factor of Brms/Bbend- Bascd on the idea that all characteristic 
scales are fixed by equating the curvature in the bending regions 
to 1/ /9,-,bend. wc havc Constructed a field-reversal-scale estimate 
of order 10 pc (Schekochihin et al. 2005). However, our ar- 
gument assumed that the field structure is exactly the same as 
in the case of isotropic viscosity with just the magnetic cutoff 
scale undetermined. In order to justify this assumption, or to 
learn otherwise, we will need a more quantitative theory. 

Thus, while it is not hard to envision how magnetic fields in 
clusters can reach observed strengths on a fairly short timescale, 
the full understanding of field structure requires more work. 
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